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Abstract 

We have studied the distribution offerees in gravitational systems through 
numerical experiments. Data were taken from an A'^-body simulation in an 
expanding universe. Before clustering, the distribution of random forces was 
represented as a Holtsmark distribution; the nearest-neighbor distribution is 
also shown as a comparison. The analytical and simulation distributions are 
in good agreement. When clustering becomes strong, the simulation result 
showed that the contribution of the forces acting on each galaxy, which is 
generated from all other galaxies, is almost entirely due to the gravitational 
attraction of its nearest neighbor. This implies that nearest-neighbor galactic 
encounters may play the main role in the dynamics of galaxy clustering. 
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1. Introduction 



One of the fundamental problems in galactic dynamics is to analyze of 
the nature of the forces acting on a galaxy that is a member of a system of 
galaxies in the large-scale structure of the universe. The average force acting 
on each galaxy in gravitational systems, such as clusters of galaxies, can be 
naturally decomposed as 



Here, is a slowly varying part, derivable from a potential function, 
which is due to the smoothed-out distribution of the galaxies in the system. 
The random force, F, is due to a relatively rapid statistical fluctuation of the 
immediate local neighborhood. As a consequence of fluctuations, the actual 
value of F at any particular moment depends on the instantaneous complex- 
ion of the local distribution of galaxies. In this sense is a stochastic force, 
which can be described through the assignment of a force-distribution func- 
tion P{F), just like any stochastic variable. The probability of a stationary 
force distribution in gravitational systems was flrst given by Chandrasekhar 
and von Neumann (1942) for an inflnite homogeneous system. It includes two 
basic assumptions. First, the stars surrounding the test star on which the 
force is being calculated are distributed with a uniform probability density. 
This implies a constant density in the local neighborhood, with no correla- 
tions among the positions of the stars. Second, it allows the number of stars 
to tend to inflnity while keeping the density constant. This force-probability 
distribution is given by the Holtsmark distribution. 

Under the above assumptions, Chandrasekhar (1943) obtained a distribu- 
tion law of the nearest neighbors in a random distribution of stars. It denotes 
the probability of force which comes from the nearest-neighbor stars, which 
is called the nearest-neighbor distribution. Chandrasekhar also showed an- 
alytically that for \F\ oo the nearest-neighbor distribution is in exact 
agreement with the Holtsmark distribution. The physical meaning of this 
agreement is that the dominant contribution to F is made by the nearest- 
neighbor galaxies. 

Attemps to examine the probability of the force distribution in an inho- 

mogeneous system were begun by Kandrup (1980, 1981). He showed that 
the distribution of random force in an inhomogeneous system is a generaliza- 
tion of the Holtsmark distribution. Antonuccio-Delogu and Atrio-Barandela 



T = 



K + F. 




2 



(1992) derived the force distribution in weakly clustered systems, where the 
two-point correlation function ^(r) is << 1. 

A numerical experiment performed by Del Popolo (1996) shows that 
the analytical results of Kandrup (1980) and Antonuccio-Delogu and Atrio- 
Barandela (1992) give a good description of the random forces in inhomoge- 
neous and weakly clustered systems. 

In the present work, we considered the evolution of the force distribu- 
tion function which comes from all the galaxies in the system (hereafter 
Total- Force Distribution), and from the nearest galaxies (hereafter Nearest- 
Neighbor Distribution) before and during gravitational clustering of galaxies 
in an expanding universe, through a numerical experiment. The motivation 
of searching for the contribution of the force acting from the nearest galax- 
ies is to examine the influence of the nearest-neighbor galactic encounters in 
gravitationally clustered systems. 

Before clustering, it was assumed that galaxies are distributed under a 
Poisson distribution. In this case, we make comparisons from both the an- 
alytical and simulation results. Analytically, we applied the Holtsmark law 
and the nearest-neighbor law. During clustering, all 1000 particles, each of 
which represents a galaxy, are examined for a Q = 1.0 universe (Peebles, 
1980) for timescales of a/ao — 4.0, 8.0, 16.0, 22.63, 32.0, where a is the scale 
factor. 

2. Total-Force Distribution and Nearest-Neighbor Distribution 

The stochastic force F was obtained by subtracting the mean force from 
the total force at the i-th galaxy. 



where m denotes the mass of the field galaxy, Vij is the position of the i- 
th galaxy relative to the j-th galaxy, r\ is the position of the i-th galaxy 
to the origin, denotes the number of neighboring galaxies, and G is the 
gravitational constant. Mtot is the total mass of the system and R is its 
radius. Since the value of at a given time depends on the instantaneous 
positions of all other galaxies, it is subject to fluctuations as these positions 
change. Even if the galactic distribution is stationary and homogeneous, 
T fluctuates around an average value due to local Poisson fluctuation in 




i?3 



(2) 
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the number density of neighboring galaxies. The fluctuating part of T that 
was indicated as F can be represented by the probability distribution of the 
stochastic force method P{F). The probability for the occurence of F in the 
range F+dF is 

P{F)dFJFydF^ = P{F)dF. (3) 

An explicit form of this distribution is shown as the Holtsmark distribution 
P{F) (Chandrasekhar, von Neumann 1942), 

P(|F|) = 4^ r eM-i.xQI\F\f'^]x&mxdx, (4) 
7r|i^| Jo 

where the normal field Q is defined by 

Q = 2.6031G'(m3/2n)2/3 (5) 

{n denotes the average number of particles per unit volume). This equation 
was applied to describe the force distribution analytically in the initial condi- 
tion when galaxies are randomly distributed. In particular, the value of |i^| 
that has the maximum probability of occurence is found to be ~ 1.6 Q. The 
asymptotic behavior of the distribution function in a strong field is found by 

P{\F\) ~ 27rG^/W/2n|F|-5/2 ^ oo). (6) 

The law of distribution of the nearest neighbor is given by 

P{r)dr = exp(-47rr^n/3)47rr^ndr. (7) 

P{r)dr denotes the probability that the nearest neighbor to a particle occurs 
between r and r -\- dr. Since in a first-neighbor approximation 

|F| = Gm/r^, (8) 

we readily obtain 

P{\F\)d\F\ = e^^[-AT:{Gmf'^nl?>\Ff^]2Ti{Gmf/MFf^d\Fl (9) 
which has an asymptotic behavior at strong fields, as 

P{\F\) ~ 27r(G'm)=^/2^|F|-5/2 ^|^| ^ (^g) 
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It can be seen to be in exact agreement with formula (6), derived from the 
Holtsmark distribution. Analytically, we can compare both the Holtsmark 
distribution and the nearest-neighbor distribution by using equations (4) and 
(9). 

A simulation of the total-force distribution can be obtained by using 



N 



N 



(nxj -y^i^ f , ' S +{Nyj - E . ' \S +{Nzi - E , . + )" 



1/2 



.ill 



Here, we use the assumption that G = m = R = 1. N is the number of 
galaxies, Fj is the stochastic force acting on the j-th galaxy, which was 
generated by all of its neighbors. The nearest-neighbor distribution was 
simulated by using 



IF-I = 



/ Xj-Xn \2 / Vj-Vn \2 / Zj - Zn V 
^ \ Xj — Xn 1^ ^ I-?/,- — 7/„ 1^/ V I — 1^/ 



I 



1/2 



(12) 



■\yj-yn\ 

where y„, Zn) is the coordinate of the nearest neighbor of the j-th galaxy. 
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3. Results 

The distributions of the forces before clustering are presented in figure 1 
and figure 2. They show the Holtsmark distribution as well as the nearest- 
neighbor distribution for both the analytical curve and the simulation result. 
The Holtsmark distribution was produced by a stronger force than that in 
the nearest-neighbor distribution. In the simulation, 1000 particles were ini- 
tially randomly distributed inside a sphere of radius R — 1. Each particle 
represented a galaxy. We can see that the simulation result gives good agree- 
ment with the analytical one. From this result we sec that we may inspect 
the distribution of forces in clustered states with 1000 particles. 



Figure 1 Figure 2 



Figure 3 shows projected distributions of galaxies during clustering for 
six timescales of a/ao, where a is the cosmological expansion parameter with 
an initial value of Oq. These were all projected views onto the X-Y plane. 
The maximum comoving radius of the simulation is always scaled to i? = 1. 



Figure 3 



Figures 4-9 present the evolution of the simulation result of both the 
total-force distribution and nearest-neighbor distribution. We used equations 
(11) and (12) for a simulation of the force distribution, which was applied 
for A'"=1000 particles. For each particular value of a/oo, we notice that both 
distributions are very similar in regions of strong force, except in regions of 
weak force. This demonstrates that if a galaxy experiences large forces, they 
are almost certainly due to the single nearest-neighbor galaxy. The region 
in which the two distributions differ most markedly is when the force is very 
weak. The frequency of occurence of the nearest-neighbor distribution within 
the weak force is higher than that of the total force distribution. 

Part b of figures 4-9 illustrates the evolution of the distribution in a 
strong-force region. As the value of a/a^ becomes bigger, the total-force dis- 
tribution becomes more similar to the nearest-neighbor distribution. After 
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several time scales of expansions, when clustering is strong, the contribu- 
tion of the force generated from the nearest neighbor of the galaxy becomes 
stronger. The force acting on the galaxies per unit mass is almost entirely 
due to the gravitational attraction of the nearest neighbors. It is also noted 
that the frequencies of occurrence of forces are nearly constant over a wide 
range of the force strength. 



Figure 4 



Figure 5 



Figure 6 



Figure 7 



Figure 8 



Figure 9 



4. Conclusions 

The result of a simulation of the force distribution before clustering is in 
good agreement with the analytical distribution. Therefore, we can use the 

data from an A^-body simulation for examining the evolution of the nearest- 
neighbor and total-force distribution, although there is as yet no analytical 
function with which to make a comparison. 
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In order to follow the influence of galactic encounters, we examine the evo- 
lution of a nearest-neighbor distribution, which is compared with the total- 
force distribution. Figures 4-9 show that as the evolution of clustering is 
in progress, the total-force distribution and the nearest-neighbor distribution 
agree over most of the range of F, especially for stronger fields. The physical 
meaning of this result is that the force acting on galaxies per unit mass is 
almost entirely due to the gravitational attraction of their nearest neighbors. 

We conclude that the influence of the nearest-neighbor galactic encoun- 
ters becomes dominant when clustering is strong. The force acting on each 
galaxy in gravitationally clustered systems in the expanding universe is al- 
most entirely due to the gravitational attraction of the nearest neighbor. 
Nearest-neighbor galactic encounters might well play the main role in the 
dynamics of galaxy clustering. If so, this may have two implications. First, 
in a gravitationally clustered system, such as clusters of galaxies, each en- 
counter could be treated as a two-body encounter representing the perturber 
galaxy and its nearest-neighbor galaxy. Second, the collisionless A^-body sim- 
ulation method (Binney, Tremaine 1987) may not be well suited for studying 
the dynamics of galaxy clustering, in which nearest-neighbor galactic encoun- 
ters play the main role, because this method assumes that the force from the 
nearest neighbor is negligible and the collective force is dominant. 

We are grateful to Dr.T.Tsuchiya for reading the manuscript and making 
useful comments. E.A. thanks all members of Department of Astronomy, 

University of Kyoto, for their hospitality during her stay. She is indebted to 
the Ministry of Education, Science, Sport and Culture of Japan for financial 
support. 
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Figure Captions 



Fig. 1. Holtsmark distribution in the analytical curve and simulation 
result. 

Fig. 2. Nearest-neighbor distribution in the analytical curve and simula- 
tion result. 

Fig. 3. Six projected distributions for time scales ofa/co = 1.0, 4.0, 8.0, 
16.0, 22.63, 32.0. 

Fig. 4. Total force distribution (solid line) and nearest-neighbor distri- 
bution (short dashes line) at a/ag — 1.0 (a) in whole; (b) enlargement of the 
large-force region. 

Fig. 5. Total force distribution (solid hue) and nearest-neighbor distri- 
bution (short dashes hue) at a/oo = 4.0 (a) in whole; (b) enlargement of the 
large- force region. 

Fig. 6. Total force distribution (solid line) and nearest-neighbor distri- 
bution (short dashes line) at a/ao = 8.0 (a) in whole; (b) enlargement of the 
large-force region. 

Fig. 7. Total force distribution (solid hue) and nearest-neighbor distri- 
bution (short dashes line) at a/ao — 16.0 (a) in whole; (b) enlargement of 
the large-force region. 

Fig. 8. Total force distribution (solid line) and nearest-neighbor distri- 
bution (short dashes line) at a/ao — 22.63 (a) in whole; (b) enlargement of 
the large-force region. 

Fig. 9. Total force distribution (solid line) and nearest-neighbor distri- 
bution (short dashes line) at a/ao — 32.0 (a) in whole; (b) enlargement of 
the large-force region. 
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